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The thickness evolution of in-plane magnetization rever- 
sal in ultrathin films is studied with a theoretical model that 
takes account of surface roughness typical of epitaxial growth. 
Guided by Neel's model, step edge sites of monolayer-height 
islands are assigned a two-fold anisotropy in addition to a 
four- fold anisotropy at all sites. Coercivity is found to de- 
pend essentially on both the film thickness and the partial 
coverage of the topmost layer. Its qualitative features are de- 
termined primarily by sample geometry and the size of the 
step anisotropy compared to the domain wall energy. Mag- 
netostatic interactions change the results quantitatively, but 
not qualitatively. Their effect is understood by comparing 
calculations with and without their inclusion. 

PACS numbers: 75.70.-i,75.60.-d 

The potential for novel physics and exciting applica- 
tions has motivated many studies of ultrathin magnetic 
films jT|. To explore new physics, it is usual to focus on 
simple model systems and equilibrium properties such as 
exchange, anisotropy, and the thermodynamic phase di- 
agram in the space of temperature and thickness. To ex- 
ploit new applications, it is typical to study more complex 
systems and non-equilibrium properties such as hystere- 
sis, domain wall motion, and magnetotransport. Com- 
mon to both is the observation that variations in surface 
roughness and film morphology often have a significant 
effect on measurement results. This observation moti- 
vated the theoretical work reported below. 

Magnetometry and the surface magneto-optic Kerr 
effect (SMOKE) || are widely used to probe magneti- 
zation reversal and the origin of coercivity in ultrathin 
films. A typical experiment reports the coercive field as 
a function of total deposited material using a thickness 
scale calibrated in small fractions of a monolayer. While 
it is generally appreciated that the films in question ex- 
hibit surface roughness, the consequences of this fact are 
not often addressed explicitly. An exception is a theoret- 
ical argument presented a few years ago by Bruno et al. 
j| Making simple assumptions regarding thickness fluc- 
tuations and the nature of domain wall pinning in films 
with perpendicular magnetization, he derived a coercive 
field He oc i~ 5 / 2 where t is the film thickness. Exper- 
imental tests of this prediction for the Co/Pd(lll) || 
and Co/Pt(lll) || systems have yielded contradictory 
results. 

In this paper, we study hysteresis and coercivity as 
a function of thickness for a model ultrathin film with 
in-plane magnetization and surface roughness typical of 




FIG. 1. Schematic view of the rough ultrathin film mor- 
phology used in this work. The indicated island geometry is 
repeated periodically. Arrows indicate local anisotropy axes 
for K 4 > and K 2 < 0. 

as-grown samples. No universal law is found. Instead, 
the coercive field is found to depend on both the film 
thickness and the partial coverage of the topmost incom- 
plete layer. Our conclusion is based on zero-temperature 
numerical simulations of magnetization reversal using a 
generalization of a classical spin model of an ultrathin 
film introduced previously 0. 

The model cubic film (Figure 1) sits atop a non- 
magnetic substrate and takes the form of a periodic ar- 
ray of square, monolayer-height, magnetic islands of side 
length L and center-to-center separation D arranged on 
top of 1 < t < 5 complete magnetic layers. To model cen- 
tered cubic lattices, the variable t is measured in units 
of to = a/2 where a is the in-plane lattice constant. 
Classical Heisenberg spins at each surface site i point 
in the direction Si. At each site in the surface, the spins 
are locked throughout the thickness of the film giving 
a two-dimensional model. This two-dimensional model 
is appropriate if the film thickness is much less than 
the exchange length. Each spin is subject to nearest- 



neighbor ferromagnetic exchange J, an out-of plane sur- 
face anisotropy K z > 0, a four-fold in-plane anisotropy 
K4 > 0, Zeeman energy from an external field H, and 
a two-fold anisotropy K2 at island perimeter sites only. 
A Neel model analysis pi shows that the latter arises 
because step edges break translational symmetry in the 
surface plane. We also take account of the fact that spins 
at step edges experience reduced exchange compared to 
spins elsewhere on an island. 
The magnetic energy is 
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where U is the film height at site i in units of to, b; 
is a unit vector parallel to the local step edge, = 
Jmin[tj, tj], and /i = [ioa 2 toMs where (Xo is the magnetic 
permeability in SI units and Ms is the saturation mag- 
netization. Note that, although K2 > corresponds to 
an anisotropy axis parallel (K2 > would be perpendic- 
ular) to the step edge, for the highly symmetrical island 
geometry studied here the energy (1) is invariant to a 
change in the sign of K2 and a simultaneous rotation of 
the external field and all spins by 90° in plane. Thus, the 
same reversal dynamics is expected for both signs of K^. 
Typical values of the material parameters, J ~ 10~ 21 
J and K4 ~ 1 x 10~ 3 mJ/m 2 (1 mJ/m 2 =l erg/cm 2 ), 
§ imply a domain wall width W ~ 8 y / J/2K 4 ~ 600a. 
The numerical results reported below all use the values 
\K 2 \ =Kz~l mJ/m 2 [|, a = 0.3 ran, and (i ~ 2x 10~ 23 
J/T in addition to those quoted above for J and K4. 
The lengths L and D are measured in spin block units of 
10a ~ W/60, a distance over which no appreciable spin 
rotation occurs. 

Beginning with a large positive value of H = Hi \\ 
[100], the local minimum of (1) was followed as the field 
was reversed adiabatically by a combination of conju- 
gate gradient minimization and relaxational spin dynam- 
ics J7J. The solid curves in Figure 2 illustrate the vari- 
ation of the computed coercive field He with coverage 
= L 2 /D 2 when the incomplete magnetic layer sits atop 
t — 1, 3, 5 complete magnetic layers. Note that He is nor- 
malized to the Stoner-Wohlfarth value Hsw = 8a 2 K4/ ' \x 
where easy-axis reversal occurs for a single domain sys- 
tem with four- fold anisotropy jl(|. Because Kz is large 
enough to keep the spins in plane, the t = l curve agrees 
with our previous XY- model results 0. They also agree 
semi-quantitatively with SMOKE data reported by Buck- 
ley et al. Q for the Cu/Co/Cu(001) system. The rather 
different t = 3 and t = 5 results can only reflect the 
decreasing importance of the surface anisotropy, which 
scales like = K^jt relative to other terms in the mag- 
netic energy. It is apparent from Figure 2 that the coer- 
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FIG. 2. Scaled coercive field as a function of L/D when 
the incomplete magnetic layer sits on t = 1,3,5 complete 
magnetic layers. D = 64. The solid curves give the re- 
sults without magnetostatics. The heavy dashed curve and 
light dashed curve give the results with magnetostatics for 
the cases of step anisotropy parallel and perpendicular to the 
steps repsectively. 

civity depends not only on film thickness but also on the 
partial coverage of the topmost incomplete layer. The 
reasons for this are best understood by study of the cor- 
responding magnetization curves. 

The hysteresis loops calculated during adiabatic cy- 
cling of the field for t — 1, 3, 5 and relatively large islands 
(G ~ |) are shown in Figure 3. The t=l reversal sce- 
nario is understood 0,|l^] ■ Reversal initiates at those step 
edges where the anisotropy axis is along the y axis, (per- 
pendicular to the saturated state spin direction) because 
spins at those sites experience the greatest torque. This 
first deviation from the positive saturated state occurs 
at H = Hn > 0, the nucleation field. The rotated spins 
form small lens-shaped domains pinned at the islands 
edges within which, near H — 0, the spins are oriented 
nearly 90° to the original saturation direction. The do- 
main walls depin from the steps at the island edges at 
a (negative) field \H\ = H s ~ (W/L)H SW || and the 
magnetization jumps to a configuration where the 90° 
state covers the islands except near the complementary 
step edges where the anisotropy axis is along the x axis, 
i.e., where the two-fold axis locally pins the spins parallel 
to the saturated state spin direction. Coherent rotation 
of the island spins occurs as the external field is made 
more negative (the coercive field occurs in this interval) 



- 1 — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 




t = 3 




-0.5 0.0 
H/H 



SW 



FIG. 3. Calculated film magnetization projected onto 
H || x for films of thickness t = 1, 3, 5 and large islands with 
L = 56 and D = 64. 
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FIG. 4. Calculated film magnetization projected onto 
H || x for films of thickness t = 1, 3, 5 and small islands with 
L — 16 and D = 64. 



until at \H\ = Ht |13|, another magnetization jump oc- 
curs when the 90° single-domain state on the island ter- 
races suffers a Stoner-Wohlfarth instability to complete 
reversal (the 180° state). A final jump occurs when the 
external field induces the complementary step-edge spins 
to flip from 0° to 180° and the reversal is complete. 

The t = 3 loop in Figure 3 differs from the t = 1 loop in 
three important ways: (1) the nucleation field is negative 
rather than positive and coincides with a jump in mag- 
netization; (2) the coercive field is smaller; and (3) there 
are two magnetization jumps rather than three. Since 
nucleation occurs at two parallel step edges per island, 
the origin of the change in sign of Hjy can be under- 
stood from the authors' recent analysis of magnetic re- 
versal on vicinal surfaces Jl4|] . We showed there that, for 
a single step (infinitely separated from other steps) on 
an otherwise flat ultrathin magnetic film, reversal nucle- 
ates at H™ = H SW (K. 2 - 1) where /C = aK^/2^/2JKl is 
a dimensionless measure of the effective step anisotropy. 
Since K, ~ 1 for realistic values of the magnetic param- 
eters (as used here), it is not surprising that H N can 
change from positive to negative as the value of de- 
creases. In fact, \Hn\ > Hs in this instance, so nucle- 
ation is accompanied by an immediate jump to the 90° 
state on the islands. 

Coherent rotation begins when the external field is 
made more negative as in the t = 1 case but the rotation is 
abruptly cut off here by a jump of the island terrace spins 
to the 180° reversed state. The corresponding jump oc- 



curred at a much larger (negative) value of external field 
for the t = 1 film. This occurs because rotation of island 
terrace spins away from the 90° state costs exchange en- 
ergy near every island edge due to the two-fold anisotropy 
at each step. The reduction in Ht seen in Figure 3 from 
t = 1 to t = 3 is thus expected because the effective two- 
fold anisotropy is reduced for thicker films. It happens 
that He = Ht in this case. Only a coherent rotation 
of the complementary step spins is required to complete 
reversal now because K% is here too small to strongly pin 
those spins in the original saturation direction. 

The t — 5 hysteresis loop in Figure 3 is easily explained 
because the formula quoted above for H^ predicts an 
even larger value for the nucleation field as K% continues 
to decrease. A simple square loop is found because |-ff/v| 
is greater than both Hs and Ht- He = Hm and the 
entire film acts as a single domain particle. 

Figure 4 is the same as Figure 3 except that the mag- 
netization curves are shown for relatively small islands 
(0 ~ 1/16). The t — 1 loop evidently differs drastically 
from the corresponding loop for large islands. The nu- 
cleation field Hn > but the much smaller island size 
(L) drives the value of Hs ~ (W/L)Hsw so large that it 
exceeds Ht- A direct jump to nearly complete reversal 
occurs at He = Hs ■ The change in sign of Hn explained 
above occurs between t=l and t — 3 as before. But now 
the sequence \Hn\ > Hs > Ht guarantees that only a 
square loop will be found for t — 3 in Figure 4. The ap- 
proach of He — Hn to Hsw observed for the t = 5 loop 



follows from the same argument based on as was 
advanced earlier for the large island case. 

Since He = Hn for t > 5, the t = 5 solid curve for He 
in Figure 2 actually shows the island size dependence of 
Hn- The results are comparable to the nucleation field 
dependence on step separation for a periodic array of 
steps separated by a finite distance Q. Those results, 
presented in Figure 7 of Ref . |l4| , show Hn rapidly ap- 
proaching zero (from negative values) as the step spacing 
decreases through values comparable to the distance be- 
tween adjacent steps used in the present work. Qualita- 
tively, nucleation is easier when steps are close together 
because the intervening spins tend to rotate also and thus 
reduce the total exchange energy cost. 

It has been argued that magnetostatics has a negligible 
effect on the hysteresis of ultrathin films with in-plane 
magnetization |15|] . We checked this explicitly by adding 
a dipole-dipole energy term to the magnetic energy (1). 
In the continuum limit, this term can be written 

E d = -yJ drM ' H D (2) 

where the dipolar field Hd is found by solving the 
Maxwell equations V • H_d = — V • M and V x Hd — 0. In 
practice, these equations are discretized and M; = MsSi 
throughout the thickness of the film at site i. In a sys- 
tem without step edges any form of magnetization inho- 
mogeneity is disfavored by magnetostatics because the 
uniform state absolutely minimizes the dipolar energy. 
In systems with step edges spin alignment parallel to 
steps is favored to avoid exposed magnetic poles. This 
breaks the symmetry noted earlier regarding the sign of 
K 2 . Roughly speaking, magnetostatics tends to enhance 
the step anisotropy when K 2 < (easy axis parallel to 
the steps) and weaken the step anisotropy when K 2 > 
(easy axis perpendicular to the steps). We examined 
both cases. 

The computational methodology used by us to eval- 
uate the dipole energy is similar to that described by 
Mansuripur |Q. The curves with short (long) dashes 
in Figure 2 represent the coercive field including mag- 
netostatics for the case of K 2 > (K 2 < 0). Practi- 
cally no effect is seen for t = 1 as suggested in Ref. [jlji) . 
But the t = 5 curves show just the combination of effects 
expected qualitatively. Nucleation at steps unavoidably 
generates magnetization inhomogeneities so a larger driv- 
ing force (Hn) is required to initiate reversal. This effect 
is partly compensated when K 2 < when the intrinsic 
step anisotropy is parallel to the step. Conversely, nucle- 
ation is even more retarded when magnetostatics works 
to further reduce K$- 

When He is determined by nucleation (all of the t = 5 
curve and t — 3 up to about 9 = 3/4), the system- 
atics of the magnetostatic results can be understood 
by noting that the source of magnetostatic energy is 
V ■ M = d x M x + d y My. At saturation, M x is a constant 
and My — 0. To lowest order, M x remains constant at 



nucleation so only the variations of M y in the y direction 
contribute to Ed- We have seen that the magnetization 
pattern at nucleation consists of lens-shaped domains 
centered on those step edges where the local anisotropy 
axis points in the y direction |L7||. At the smallest cov- 
erages when the islands are very small, M y = M y (x,y) 
is a function of both x and y so a finite magnetostatic 
effect is expected. This effect increases as the island size 
increases initially because more step edge contributes to 
nucleation. But for large enough island size, the magne- 
tization pattern at nucleation more nearly resembles the 
vicinal surface case |T^] and My = M v (y) when K 2 < 
and M y = M y (x) when K 2 > 0. The effect of magneto- 
statics thus increases in the former case and disappears 
in the latter case as the islands grow larger. This is the 
trend seen in Figure 2 except for the very highest cover- 
ages of t = 3 when He is not determined by nucleation 

IS 

Note also that magnetostatics has a weaker effect for 
t = 5 than for t = 3. This is so because the thicker 
film experiences a relatively weaker step anisotropy and 
behaves like a Stoner-Wohlfarth particle. The magneti- 
zation is more uniform and V • M and Ed arc reduced 
accordingly. 

In summary, we have used a classical spin model to 
study zero temperature, in-plane, magnetization reversal 
in ultrathin films with surface roughness characteristic of 
the epitaxial growth process. The model includes four- 
fold anisotropy at all sites, two-fold anisotropy at step 
edge sites, and magnetostatic interactions. The coercive 
field was found to vary non-monotonically (but explica- 
bly) as a function of both film thickness and the partial 
coverage of the surface layer. We conclude that measure- 
ments of the thickness dependence of the coercive field 
and related hysteretic properties of ultrathin films with 
in-plane magnetization should not be expected to exhibit 
universal behavior. The results can depend sensitively 
on the film morphology which, in turn, can depend sen- 
sitively on growth conditions. Comparisons of nominally 
identical measurements on nominally identical systems 
in ignorance of the film morphology is unwarranted and 
can be seriously misleading. 
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